Incompatibility of modulated checkerboard patterns with the neutron scattering 

resonance peak in cuprate superconductors 
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Checkerboard patterns have been proposed to explain the real space structure observed in STM 
experiments on BSCCO and Na-CCOC. However, simple checkerboard patterns have low energy 
incommsensurate (IC) spin peaks rotated 45 degrees from the direction of the charge IC peaks, 
contrary to what is seen in neutron scattering. Here, we study modulated checkerboard patterns 
which can resolve the low frequency inconsistency. Using spin wave theory, we explore the spin 
response of these superstructures and find that the high energy response is inconsistent with neutron 
scattering results. In particular, the modulated checkerboard structures are incapable of supporting 
the experimentally well-established resonance peak at {tVjTv). 
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PACS numbers: 74.25.Ha, 74.72.-h, 75.30.Ds, 76.50.+g 



I. INTRODUCTION 



Spin and charge ordering have been topics of great in- 
terest in strongly correlated electronic systems. A key 
aspect of many strongly correlated models is that differ- 
ent terms in the Hamiltonian compete, which introduces 
a type of electronic frustration since solutions cannot be 
found which simultaneously minimize all terms in the 
Hamiltonian. These competing interactions can lead to 
spontaneous nanoscale electronic structure. Indeed, sev- 
eral locally inhomogenous electronic phases have been 
proposed, involving charge order, spin order, and orbital 
order among others, in strongly correlated materials such 
as cuprate superconductors, nickelates, manganites, and 
related perovskites. Charge order is amenable to de- 
tection through probes that directly measure charge de- 
grees of freedom, such as scanning tunneling microscopy 
(STM),—^ while the presence of spin order can be di- 
rectly detected through neutron scatteringi^i^i^iii^ Un- 
fortunately, the charge patterns which most naturally ex- 
plain the STM data have often been incompatible with 
the spin patterns which most simply explain the neu- 
tron scattering data. Part of the challenge has been that 
materials which are most amenable to STM studies (i.e. 
strongly layered materials) are least amenable to neu- 
tron scattering (which requires large crystals), and vice 
versa. Further complicating a clear and consistent inter- 
pretation of the data set as a whole is that while neutron 
scattering is a bulk probe, STM is confined to the surface. 

Recent experimental advances have made possible the 
detection of high energy neutron scattering spectra^ i^'-'i^^ . 
This has led to the discovery that the magnetic excita- 
tions in several cuprates, i.e., La2_xBa2;Cu04 (LBCO), 
La2-xSr^Cu04 (LSCO), and YBaaCuaOg+s (YBCO), 
exhibit universal behavior— One prominent fea- 
ture is that at intermediate energies, there is a resonance 
peak^** at (tt, tt) formed by the merging of the low energy 
incommensurate response with a high energy response 
whose incommensurate structure is rotated 45° from that 
of the low energy excitations. The resonance peak is ob- 
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FIG. 1: (Color online) Checkerboard pattern, a) Simple 
checkerboard with spacing p = 4. b) NoncoUinear checker- 
board with p = 4. c) Modulated checkerboard pattern with 
spacing p = 2. d) Modulated checkerboard pattern with spac- 
ing p = 3. The dotted lines represent unit cells. 



served in the pseudogap regime {i.e. the low temperature 
normal state of the underdoped cuprates), but sharp- 
ens as temperature is lowered into the superconducting 
state. The relation between the resonance peak and the 
emergence of superconductivity is still under researchii^ 
The magnetic excitations in these materials have been 
explained using stripes, a unidirectional modulation of 
spin and charge . ^"'^^i^'^^'^^i-'^^'^" . In these models, the res- 
onance peak is a saddlepoint in the dispersion. 
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FIG. 2: (Color online) Spin (black) and charge (red) expected 
diffraction peaks of checkerboards in reciprocal lattice space, 
corresponding to the real space patterns of Fig. [1] (a) , (b) , (c) 
and (d). 



STM can directly detect charge order at the sur- 
face. Checkerboard patterns (a two-dimensional mod- 
ulation of charge) have been proposed to explain the 
real space structure observed in STM experiments on 
BiaSrsCaCuaOs+^i (BSCCO) and Ca2-a;Na:,Cu02Cl22. 
(Na-CCOC). The charge modulations were character- 
ized by checkerboad patterns with spacing approximately 
4a X 4a, where a is the lattice spacing. While this 
spacing is in agreement with the incommensurability 
observed in neutron scattering, the presence of a true 
two-dimensional modulation such as a checkerboard has 
not been confirmed by neutron scattering or other bulk 
probes. 

Within microscopic models, either stripes or checker- 
boards can be stabilized by tuning parameters, such 
as the next-nearest neighbor hopping t' ^ The most 
natural checkerboard patterns to arise out of micro- 
scopic models are "simple" checkerboards, such as the 
one shown in Fig. 1(a) Such simple checkerboards are 
in fact ruled out experimentally, because the low en- 
ergy charge peaks are rotated 45° from the low en- 
ergy spin peaks, contrary to what is observed in neu- 
tron scattering Later, modulated checkerboards 
were proposed^^i^i as shown in Fig. 1(c) and 1(d) Recent 



experimental work has ruled out the modulated checker- 
boards in Lai.48Ndo.4Sro.i2Cu04, based on the relative 
intensities of the incommensurate (IC) spin peaks in dif- 
ferent magnetic Brillouin zones.— The authors of Ref. d 
also proposed a new noncoUinear checkerboard pattern 



as shown in Fig. 1(b) which is consistent with all of 



checkerboards there is no possibility of a resonance peak 
at (tt, tt), which rules out these structures in all materials 
where a resonance peak has been observed. We further 
argue that the newly proposed noncoUinear checkerboard 
also lacks a resonance peak. 



II. MODEL AND METHOD 

In this paper, we study the magnetic excitations of var- 
ious checkerboard patterns. Simple checkerboards of the 
type shown in Fig. 1(a) have been studied previously by 
us and other s — i ^ i^ 



These types of patterns, in which 
the sign of the Neel vector changes across each charge 
line (whether vertical or horizontal), are always found to 
have incommensurate (IC) spin peaks which are rotated 
45° from the IC charge peaks, contrary to what is ob- 
served in experiment from STM (which can measure the 
charge peaks), and neutron scattering (which can mea- 
sure the spin peaks, and sometimes also the charge peaks 
as well). Modulated checkerboards, like those shown in 
Figs. 1(c) and 1(d) offer a consistent description of the 



low energy data, as does the noncoUinear checkerboard 
shown in Fig. 



1(b) 



Here, we extend our previous work 
on the magnetic excitations of stripes and simple checker- 
boards to the modulated and noncoUinear checkerboard 
patterns. 

We study the spin excitations within the Heisenberg 
model, 



(1) 



where the indices i and j run over all spin sites, and 
Jij represents the spin coupling. We have assumed that 
the charge degrees of freedom can be integrated out to 
produce the effective spin couplings of the model. The 
main effect of the charge degrees of freedom is to form 
antiphase domain walls across which the Neel vector of 
the spins changes sign. 
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As shown in Fig 
neighbor couplings are Ja (antiferromagnetic 
nearest neigbor couplings across a domain wall are Jt 
(also antiferromagnetic) . We make the physically reason- 
able assumption that Ji, is small compared to Ja- The 
dotted lines in Fig. [T] show the unit cells. Note that in 
Figs. 1(c) and 1(d) the charge domain walls run diago- 
nally. In Fig. 1(c) the diagonal spacing between domain 
walls is p = 2 in units of the diagonal spacing V2a. This 
configuration has 10 spins in the unit cell. In Fig. 
the diagonal spacing between domain walls is p = 
the same units, and there are 24 spins in the unit cell. 



1(d) 



For Figs. 1(c) and 1(d) the spin ground states shown are 



the low energy data. We show here that for modulated 



unfrustrated, there is long range magnetic order, and the 
elementary excitations can be captured by the spin wave 
treatment below. Note that the noncoUinear checker- 
board is not a ground state of this spin model. We will 
return to the magnetic excitations of this state later. 

In order to study the magnetic excitations, we use the 
spin wave method to calculate the magnon excitation 
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FIG. 3: Spin wave dispersion and intensities for a modulated 
checkerboard pattern with spacing p — 2. The upper panel 
is along the (k^, it) direction and lower panel is along the 
diagonal direction {kx, kx). The energy E is in units of JaS. 



spectrum and the zero-temperature dynamic structute 
factor 



(/|5^(k)|0)p<5(c.-c.^), (2) 



where |0) is the magnon vacuum state and |/) denotes 
the final state of the spin system with excitation energy 
u)f. 5'(k, Lo) is proportional to the expected neutron scat- 
tering intensity. We show two different sizes of spins in 
Figs. 1(c) and |l(d)" We take the z component of the 
larger spins to be twice that of the smaller spins, in or- 
der to take into account that the spin modulation is not 
a square wave. However, our major conclusions arc inde- 
pendent of the details of this spin ratio. 



III. RESULTS 

We first discuss the zero-frequency response, shown in 
Fig. [51 The greyscale image shows the expected diffrac- 
tion peaks and relative intensities from the spin texture 
for each panel from Fig. [TJ The diffraction peaks from the 
corresponding charge modulation are denoted schemati- 
cally by the red circles, for fiduciary points around (0,0). 
For simple checkerboards, spin diffraction peaks are ro- 
tated 45 degrees from the direction of charge diffraction 
peaks, unlike what is seen in experiment. For the non- 
coUinear checkerboard of Fig 



1(b) and the modulated 
checkerboards of Figs. 1(c) and 1(d) the relative orienta- 
tion of the main spin and charge IC peaks are consistent 
with experiment. The noncoUinear checkerboard gives 
more satellite peaks than the modulated checkerboards, 
even when the spacing of the noncoUinear checkerboard 
is comparable to that of the modulated checkerboards, 
which may help distinguish these patterns. However, the 
higher harmonic spin peaks get weaker with increasing p. 
For the experimentally relevant case of p = 4, the main 



IC peaks are nine times stronger than the next nonzero 
harmonic peaks, which may make them difficult to de- 
tect. Note that for modulated checkerboards of even 
spacing [p — even), there is a charge diffraction peak 
located at (tt, tt), contrary to what is observed in experi- 
ment. The strength of this peak decreases as the charge 
profile is made smoother. 

Spin wave dispersions with intensities for the p — 2 
modulated checkerboard case are shown in Fig. [3] at 
Jb = O.lJa and Jb = 0.5Jq along the {k^, tt) and [k^, kx) 
directions. Note that although zero-frequency weight is 
forbidden at (tt, tt) because of the presence of antiphase 
domain walls in the Neel vector, nevertheless, the point 
(tt, tt) is a reciprocal lattice vector, and so the spin-wave 
dispersion must approach w — > O at this wavevector. 
There are several defining characteristics of the resonance 
peak, but the most well-established is that the mode oc- 
curs at finite frequency. The reciprocal lattice structure 
of the modulated checkerboard patterns therefore forbids 
the appearance of a resonance peak in the acoustic band. 

For comparison, the simple checkerboard patterns we 
studied in Ref. 20 are also incapable of supporting a res- 
onance peak in the acoustic band, but for a different 
reason. In these cases, (tt, tt) is not a reciprocal lattice 
vector, and so the acoustic mode at (tt, tt) has finite fre- 
quency. However, in the case of simple checkerboards, 
the acoustic band reaches a local maximum at (7r,7r), 
rather than the saddlepoint found in stripe phases. The 
saddlepoint structure has been shown to capture the phe- 
nomenology of the resonance peak^i^iS including the 
finite- frequency peak in the integrated spin structure fac- 
tor S{lj)^ as well as the presence of incommensurate scat- 
tering which smoothly connects to the resonance peak 
both below and above it in frequency. By contrast, while 
the simple checkerboards have a peak in S{lj) at finite fre- 
quency, there is no weight immediately above the (tt, tt) 
point in frequency. This is one of a few reasons why 
simple checkerboard patterns have been ruled out. 

For the modulated checkerboards studied here, one 
may consider the possibility of the resonance peak ap- 
pearing in one of the optical bands. There is a gap be- 
tween the acoustic band and optical bands, similar to the 
simple checkerboard configurations. The gap is sizable 
when Jb is small. The acoustic band begins to touch the 
optical bands at Jb = J a- For the physically reasonable 
assumption that Jb is smaller than J a, any weight ap- 
pearing in an optical band at (tt, tt) is too far removed 
from the incommensurate scattering at low frequency to 
be a candidate for the resonance peak. Furthermore, the 
lack of a saddlepoint structure in the optical bands fur- 
ther rules out a resonance peak-like structure (see Fig. [31 
where instead of a saddlepoint, two bands cross at (tt, tt)). 

Constant energy cut plots are experimentally measur- 
able and useful for analyzing the microscopic structures. 
Fig. [4] shows representative constant energy cut plots 
for the p = 2 structure of Fig. 1(c) 
•h/Ja = 0.1. For the modulated c 



at coupling ratio 
leckerboards, it is 



clear that the direction of the low energy IC spin peaks 
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FIG. 4: (Color online) Constant energy cuts with windows 
O.lJaS' for a modulated checkerboard pattern with diagonal 
spacing p = 2 as described in the text. The coupling ratio is 
Jb/Ja = 0.1. The energy E is in units of JaS. In each plot, 
we have integrated over an energy window AE = ±0.1Ja. 
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FIG. 5: (Color online) Constant energy cuts with windows 
0.1,/a5' for a modulated checkerboard pattern with diagonal 
spacing p = 2 as described in the text. The coupling ratio is 
Jb/Ja = 0.5. The energy E is in units of JaS. In each plot, 
we have integrated over an energy window A_E = zbO.lJa. 



are consistent with that observed experimentally for the 
low energy IC charge peaks, whereas this was not the 
case for simple checkerboards. We find that although a 
spin-wave cone must be present at low energy, the in- 
tensity of the spin structure factor is not uniform on the 
cone, as shown in the E — 0.2 panel. As in our previous 
studies of arrays of antiphase domain walls (whether in 
stripe or simple checkerboard patterns), for small cou- 
pling ratio Jb/Ja, the intensity is strongest on the inner 
branch of the spin wave cone, i.e. the side closest to 
the (tTjTt) point. We take this to be a generic feature of 
spins which are weakly coupled across arrays of antiphase 
domain walls. 

Note also the presence of a faint spin wave cone ema- 
nating from the (tTjTt) point. This spin wave cone is re- 
quired by symmetry, since the (tt, tt) point is a reciprocal 
lattice vector of the modulated checkerboard structures. 
However, due to the antiphase domain walls, weight is 
forbidden at zero frequency at (vr, tt) and this makes the 
central cone quite faint at finite frequency compared to 
those emanating from the main IC points. At interme- 
diate energies {E = 0.4), the spin wave cones touch each 
other. At higher energies, just below the top of the acous- 
tic band, incommensurate peaks are once more observed, 
now rotated to the diagonal direction, as shown in the 
E — 0.46 panel. Note the complete absence of weight at 
the (7r,7r) point, precluding a resonance peak from this 
type of modulated checkerboard. 

For comparison, in Fig. [5] we show similar constant en- 
ergy cuts, with a stronger coupling ratio Jb/Ja = 0.5. 
These have a steeper dispersion (as is evident from 
Fig. [3]), and the spin wave cones appear simply as incom- 
mensurate peaks in the lowest energy panel, E = O.A. As 
energy increases, the spin wave cones become visible as 
shown for E = 1. As with the weaker coupling case, the 
intensities are not uniform on the spin wave cone, and 



the intensities peak on the side facing {tt,7t). Note that 
at this energy, the faint spin-wave cone emanating from 
(tt, tt) becomes visible. At higher energies {E — 1.4), the 
spin wave cones merge. At yet higher energies {E = 1.6), 
the top of the acoustic band produces incommensurate 
peak structure at finite frequency, with peaks rotated 45° 
from the low energy IC peaks. As with the weaker cou- 
pling ratio, the acoustic band is forbidden by symmetry 
to support any weight at the (tt, tt) point, ruling out these 
structures as being able to support a resonance peak. 

We now briefly comment on the noncollinear checker- 
board proposed in Ref. @ and reproduced schematically 



in our Fig. 1(b) Like the modulated checkerboard pat- 



terns, this pattern has IC spin peaks which are in the 
same direction as the main IC charge peaks, consistent 
with low energy experimental data. While modulated 
checkerboards have been challenged based on the exper- 
imental results of Ref. 's', the authors were unable to 
definitively rule out the noncollinear checkerboard ( "two- 
q structure" ) . We are not able to calculate the expected 
magnetic excitation spectrum within the framework of 
the current model, since this pattern is not a valid ground 
state of the Hamiltonian we consider here (see Eqn. [1]) 
This does not preclude it being the ground state of some 
other model. What we can say about this structure is 
that to the extent that it supports Goldstone modes, it 
will have the same limitations of the modulated checker- 
board patterns discussed here. This is because the (tt, tt) 
point is a reciprocal lattice vector of the noncollinear 
checkerboard, and so it is constrained by symmetry to 
have a spin wave cone emanating from the (tTjTt) point. 
This point is also forbidden by symmetry to have any 
weight at zero frequency because the pattern has no net 
Neel vector at {tt,tt). Since the intensity of a Goldstone 
mode must be continuous in frequency, the spin wave 
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cone emanating from (tt, tt) is constrained to be quite 
weak. Therefore, like the modulated checkerboards stud- 
ied here, the noncoUinear checkerboard will be unable 
to support a resonance peak at finite frequency, except 
perhaps at unphysically high frequencies in an optical 
branch. 



IV. CONCLUSIONS 

In conclusion, we have shown that while modulated 

checkerboard patterns have low energy incommensurate 
charge and spin peaks that are consistent with STM mea- 
surements in BSCCO and Na-CCOC and with neutron 
scattering measurements in and neutron scattering ex- 
periments on lanthanum cuprates and YBCO, the finite 
frequency magnetic excitations of these structures are 
incompatible with experimental findings. In particular, 
modulated checkerboards are forbidden by symmetry to 
have weight at (tt, tt) in the acoustic branch, precluding 
the possibility of a resonance peak in this branch. Al- 



though optical modes arc not forbidden to have weight 
at (tt, tt), the structure of the optical modes around (tt, tt) 
is incompatible with the phenomenology of the resonance 
peak. We argue that similar physics constrains the acous- 
tic branch of the recently proposed noncollinear checker- 
board, relegating any possibility of a resonance peak to 
unphysically high energies. We conclude that to date 
no checkerboard pattern has been proposed which satis- 
fies both the low energy constraints and the high energy 
constraints imposed by the current body of experimental 
data in cuprate superconductors. 
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